•  •  ’  (ZD 

AD-A207  323 


OFFICE  OF  NAVAL  RESEARCH 
Contract  N00014-82-K-0576 

Technical  Report  No.  41 

ADSORPTION  AND  COADSORPTION  OF  CO  AND  NO  ON  THE  RH(100)  SURFACE: 

A  THEORETICAL  ANALYSIS 

by 

D.  L j .  Vuckovid,  S.  A.  Jansen  and  R.  Hoffma  i 

Department  of  Chemistry 
Cornell  University 
Baker  Laboratory 
Ithaca,  NY  14853-1301 

April  1989 


Reproduction  in  whole  or  in  part  is  permitted 
for  any  purpose  of  the  United  States  Government 


This  document  has  been  approved  for  public  release 
and  sale;  its  distribution  is  unlimited 


89  3  051 


BHiPiaiBaHigBP 


ilFICATION  OF  THIS  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 

unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION  /DOWNGRADING  SCHEDULE 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

#41 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Department  of  Chemistry 


6c  ADDRESS  (City,  State,  and  ZIP  Coda)  _ 

Cornell  University 
Baker  Laboratory 


iMir-w- 


uiraaitiu 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


8c.  ADDRESS  (City,  State,  and  ZIP  Code) 


11.  TITLE  (Include  Security  Classification) 


6b.  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 
(If  applicable) 

ONR 


7b  ADDRESS  (City,  State,  and  ZIP  Code) 

800  Quincy  Street,  Arlington 


8b.  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

Report  #41 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  I  PROJECT 
ELEMENT  NO.  I  NO. 


Adsorption  and  Coadsorption  of  CO  and  NO  on  the  Rh( 100) Surface:  A  Theoretical  Analysis 


12.  PERSONAL  AUTHOR(S) 

D.  Lj.  Vuckovic,  S.  A.  Jansen  and  R.  Hoffmann 


13a.  TYPE  OF  REPORT  Il3b.  TIME  COVERED 

Technical  Report  #41  |  from _ to 


14.  DATE  OF  REPORT  (Year,  Month,  Day)  IlS.  PAGE  COUNT 

April  14,  1989  | 


19.  ABSTRACT  (Continue  on  reverse  if 

adsorption  of  <&,  Ifl&and  their  coadsorption  on  the  Rh(100)  surface  ard  dYscussS? 
in  terms  of  extended  Huckel /tight  binding  calculations.  Experimentally,  it  is  not  easy 
to  resolve  the  linear  and  bent  forms  of  NO  and  for  this  reason  our  calculations  serve  to 
examine  not  only  the  proposed  adsorption  geometries  but  also  as  means  to  investigate  the/ 
NO  dissociation  mechanism.  This  dissociation  is  believed  to  be  the  fundamental  step  in/ 
the  C0/N0  reaction.  The  adsorption  site  and  geometry  are  discussed  for  several  NO  / 
chemisorptions.  In  addition,  a  discussion  of  the  C0/N0  reaction  for  the  proposed  adsorpt 
of  NO  is  presented..  From  this  analysis,  it  is  clear  that  predissociation  of  NQijnto  wit 
components  can  lead  to  recombination  producing  the  desired  reaction  products, (CO^andQtd. 


into  a 

iv>and( 


Furthermore,  the  coadsorbate  interactions  can  lead  directly  to  the  desired  products. 


20  .DISTRIBUTION/  AVAILABILITY  OF  ABSTRACT 
S  UNCLASSIFIEO/UNLIMITED  □  SAME  AS  RPT.  Q  OTIC  USERS 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Id  Hoffmann 


00  FORM  1 473, 84  MAR  83  APR  edition  may  be  used  until  exhausted. 

All  other  editions  are  obsolete. 


121.  ABSTRACT  SECURITY  CLASSIFICATION 

unclassified 


Area  Code)  22c.  OFFICE  SYMBOL 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


Adsorption  and  Coadsorption  of 
CO  and  NO  on  the  Rh(lOO)  Surface: 
A  Theoretical  Analysis 


Dragan  LJ .  Vu£kovi4* ,  Susan  A .  Jansen  and  Roald  Hoffmann* 
Department  of  Chemistry  and  Materials  Science  Center 
Cornell  Unversity,  Ithaca,  HY  14853-1301 


Abstract : 
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the  linear  and  bent  forms  of  NO  and  for  this  reason  our 

calculations  serve  to  examine  not  only  the  proposed  adsorption 
geometries  but  also  as  means  to  investigate  the  NO  dissociation 
mechanism.  This  dissociation  is  believed  to  be  the  fundamental 
step  m  the  CO/NO  reaction.  The  adsorption  site  and  geometry 
are  discussed  for  several  NO  chemisorptions .  In  addition,  a 
discussion  of  the  CO/NO  reaction  for  the  proposed  adsorption  of  NO 
is  presented.  Prom  this  analysis,  it  is  clear  that  predlssocxatlon 
of  NO  into  atomic  components  can  lead  to  recombination  producing 
the  desired  reaction  products,  CO 2  and  H£>  Furthermore,  the 
coadsorbate  interactions  can  lead  directly  to  the  desired  products. 
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The  catalytic  reduction  of  nitric  oxide  with  carbon 
monoxide  over  transition  metal  surfaces,  forming  CO 2  and  N2  as 
reaction  products ,  Is  important  in  pollution  abatement .  Among 
the  transition  metals,  rhodium  is  the  most  effective  catalyst  for 
this  purpose ,  and  in  the  automotive  exhaust  environment  it  is  the 
only  catalyst  which  works  effectively  to  reduce  HO.  Details  of 

the  catalytic  mechanism  are  still  unclear .  In  this  work ,  we  will 
study  theoretically  the  coadsorption  of  CO  and  NO  on  the  Rh  ( 1 00 ) 

surface ,  as  well  as  the  behavior  of  each  individual  adsorbate . 

Let  us  first  review  briefly  some  previous  experimental 
and  theoretical  work  In  this  area .  The  adsorption  of  CO  on 
transition  metal  surfaces  has  been  studied  extensively1*10.  The 
experimental  work  has  been  accompanied  by  much  theoretical 
consideration11'10.  On  early  transition  metals  CO  chemisorbs 

dlssoclatlvely ,  and  on  late  transition  metals ,  molecularly 1 » 1 2 . 
On  most  rhodium  surfaces,  CO  is  adsorbed  molecularly  with  the 
carbon  end  down .  The  adsorption  geometries  and  site  preferences 
for  low  index  Rh  surfaces  are  well  known3'0.  At  lower  coverages, 
the  only  CO  state  present  is  a  molecular  species  in  an  on -top 
site,  i.e.  each  CO  being  bound  to  a  single  Rh  surface  atom. 

With  Increasing  CO  coverage,  a  new  compressed  surface  pattern 
is  formed,  with  some  CO  starting  to  occupy  two -fold  bridging 
sites3'0*17.  The  simplest  general  picture  of  CO  chemisorption 
postulates  electron  donation  from  the  CO  5c  (HOMO)  orbital 
into  suitable  transition  metal  d  orbitals  and  back  donation 
from  the  surface  d  orbitals  into  CO  2*"  (LUHO)  orbitals. 

This  is  usually  referred  to  as  the  Bly holder  model10.  There 
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is  some  evidence  that  this  simple  picture  may  not  he  complete . 
Other  CO  orbitals  may  be  important ,  especially  the  4o 1 6  * 1 9 . 
For  a  previous  theoretical  study  of  the  CO  chemisorption  on 
transition  metal  surf  aces  in  this  group  see  ref .  l  a . 

Nitric  oxide  chemisorption  on  transition  metal  surfaces 
has  also  been  studied  extensively 1  »*•  *20“2 5,  it  is  much  more 

complex  and  not  as  well  understood  as  CO  chemisorption.  The  NO 

ligand  in  both  discrete  complexes  and  on  surfaces  appears  to  have 
the  freedom  to  bond  linearly  or  bent.  In  nitrosyl  transition 
metal  complexes.  the  extremes  can  be  represented  as  H-(NO)*,  or 
H-  (HO )  ” .  with  the  first  leading  to  a  linearly  bonded  and  the 
latter  to  a  bent  NO  configuration  (NO  axis  being  substantially 
tilted  away  from  the  surface  normal ) 25-29 .  in  general,  NO 

undergoes  dissociative  chemisorption  more  readily  than  CO.  For 
NO  on  group  VIII  metals  there  Is  evidence  for  low  temperature 

dissociative  adsorption,  as  well  as  molecular  adsorption  in  a 
variety  of  geometries,  including  a  bent  form,  a  bridged  form  and 
a  linear  form30.  However,  unambiguous  evidence  for  bent  forms  is 
difficult  to  obtain  because  of  the  complexity  of  the  NO  spectra. 

The  adsorption  and  decomposition  of  nitric  oxide  on 
rhodium  has  received  some  attention  in  the 
literature2  5 » 20 » 23 » 25 » 31  »32.  Both  molecular  and  dissociative 
chemisorption  is  found ,  with  the  latter  leading  eventually  to  Og 
and  Hg  desorption.  Just  to  illustrate  the  complexity  of  NO 
adsorption  on  Rh  surf  aces ,  it  is  worthwhile  to  mention  that 
Dubois  et  al .  5  have  reported  at  least  five  different  NO  species 
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on  Rn  ( 1 1 1 )  using  man  resolution  EELS .  For  a  detailed 
theoretical  study  of  HO  on  Hl(iii),  see  another  contribution  from 
this  group33* 

The  CO/NO  reaction  on  the  transition  metals  has  been 
investigated,  for  instance,  on  clean  poly  crystalline  Pt34,  Pt(ill)35, 

Pt ( 1 00 ) 5  and  Ru ( 001 ) 35 .  in  general,  it  is  accepted  that  the  NO 
dissociation  is  a  prerequisite  for  the  reaction  to  occur,  at 
least  for  some  surf  aces .  On  surfaces  such  as  Pt  ( 1 1 1 ) ,  in  which  HO 
is  adsorbed  only  molecularly,  no  reaction  with  CO  occurs34*35. 

The  mechanism  of  nitric  oxide  reduction  is  unclear 
and  despite  its  significance  there  have  been  relatively  few 
mechanistic  studies  examining  the  CO/HO  process  on  well  defined 
rhodium  surfaces.  The  few  that  exist  are  recent:  Rh(ilO)4,  Rh(33i)5, 
polycrystalline  Rh 37 ,  Rh  ( 1 1 1 ) 35  and  Rh  ( 1 00 ) 39 .  There  have 
been  numerous  investigations  of  the  CO/NO  reaction  over 
supported  rhodium  such  as  Rh/AlgOj20*40  and  Rh/Slo241 .  The 
basic  question  about  the  reaction  mechanism  is  whether  the  major 
steps  involve  HO  decomposlton ,  or  if  it  is  a  true  bimolecular 
reaction  between  adsorbed  CO  and  HO .  For  the  Rh  ( i  i  l )  surface ,  it 
wasrecently  shown 3  ^  that  the  mechanism  goes  exclusively  through 
the  NO  decomposition  step. 

In  an  attempt  to  shed  some  light  on  NO  adsorption  and  its 
coadsorption  and  reaction  with  CO  we  present  here  the  theoretical 
investigation  of  the  CO/HO  system  on  a  Rh  ( 1 00 )  surface ,  using  the 
extended  Htickel/tight  binding  method42.  This  procedure  has  well- 
known  deficiencies  but  in  general  describes  simply  the  reactive 

trends  responsible  for  surface  chemistry . 
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Clean  Rh(lOO) 

The  Rh(100)  single  crystal  surface  is  well  defined  and 

nonreconstructed2 » 31 » 39 *43,  giving  a  sharp  (ixl)  LEBD  pattern. 
For  this  reason ,  the  model  chosen  for  the  calculation  was  a  three 
layer  two  dimensional  slab  of  the  fee  rhodium  metal  structure44 
with  lattice  constant  of  3.d032A  and  nearest  neighbor  separations 
of  2.0093A  (1).  Whether  the  use  of  the  three  layer  slab  model 
to  represent  the  surface  in  our  calculation  is  appropriate  has , 
been  discussed  earlier  12,33,45,40. 


top 


1 


surface 

inner 

surface 


side 


The  calculated  density  of  states  for  a  three  layer  slab 

shows  a  d  band  between  -7.0  and  -13.0  eV ,  and  dispersed  s  and  p 
bands  starting  from  -12.0  eV  and  up  to  20  eV .  The  Fermi  level  is 
at  0.55  eV  (Figure  la).  Since  the  surface  atoms  have  a  lower 
coordination  number  than  the  bulK  atoms,  their  states  are  less 
dispersed  than  the  bulk  states.  If  the  Fermi  level  falls  above 
the  midpoint  of  the  d  block ,  as  it  does  for  rhodium ,  then  the 
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surface  should  have  a  negative  charge  with  respect  to  the  bulk46. 

In  our  case,  the  total  charge  on  the  surface  rhodium  atom  is  an 
exaggerated  9.32  electrons  and  the  bulk  atom  have  a  total  of 
6 . 35  electrons . 

CO  on  the  RhUOOl  Surface 

CO  readily  chemisorbs  on  the  Rh(l00)  surface 
with  a  sticking  coefficient  near  unity .  The  chemisorption  is 

molecular  and  no  evidence  for  dissociation  was  observed ,  as  for 

all  late  transition  metal  surfaces1.  On  the  Rh(lOO)  surface,  CO 
site  preference  is  coverage  and  temperature  dependent2 » 47 .  It 
is  necessary  to  adsorb  CO  gas  at  25°C  or  below  to  obtain  well- 
ordered  surface  structures ,  with  disordering  temperature  about 
125°C.  At  lower  exposures,  a  c(2x2)  adsorption  LBED  pattern  was 
observed  and  at  higher  exposures  above  1  L  this  structure  is 
compressed  into  a  hexagonal  overlayer ,  which  is  a  split  ( 2x  l ) 
pattern .  The  saturation  coverage  is  about  e  =  0 . 62 .  At  lower 

temperatures  (-100  K)  some  other  complicated  structures  were 
observed47  similar  to  the  coincidence  lattice  pattern  reported  by 
Tucker40,  with  a  rectangular  unit  cell  which  is  4xl  times  the  Rh 
face -centered  cell.  This  result  was  later  questioned  by  Castner 
et  al . 2  since  it  was  unclear  whether  these  studies  were  carried 
out  on  an  initially  clean  surface .  Kim  et  al . 47a  have  found  a 
similar  compressed  structure. 

In  general ,  we  can  say  that  at  lower  exposure  a  regular 
c (2x2)  pattern  (coverage  of  0.5)  is  observed  with  the  on -top 
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*  adsorption  geometry.  Higner  exposures  lead  to  a  compressed  leed 
structure  with  a  saturation  coverage  near  0.0,  involving  a  new 
adsorption  state  with  a  Higher  coordination  number.  This  can  be 
assigned  to  a  two -fold  bridging  state39.  The  difference  of  the 
neats  of  desorption  between  these  two  states  was  calculated  by 

Hendershot  and  Hansen39  to  be  3.2  Heal /mol. 

The  molecular  orbitals  of  an  isolated  CO  are 
well  Known49.  The  highest  occupied  molecular  orbital  (HOHO)  is 
mainly  the  carbon  lone  pair  5 c  ( 2a ) ,  and  its  energy  is  lower  than 
the  d  states  of  most  metals.  In  the  generally  accepted  picture  of 
CO  chemisorption,  5a  is  stabilized  through  interaction  with 
surface  d  states  relative  to  other  CO  levels.  The  lowest 
unoccupied  molecular  orbital  (LUHO)  of  the  free  CO  molecule 
the  2s"  consists  of  two  antibonding  combinations  (21>)  with 
larger  coefficients  on  the  carbon  atom.  As  we  have  already 


METAL 

METAL 

METAL 

a 

2 

b 

mentioned,  these  two  orbitals  are  mainly  involved  m  the 
chemisorption  mechanism  (Bly holder  model ) .  Upon  chemisorption 
the  2v”  derived  orbitals  are  partially  filled  and  consequently , 
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the  CO  bond  is  weakened.  Aside  from  the  frontier  orbitals, 
other  CO  orbitals  are  usually  not  taken  into  account  because 
they  are  too  distant  in  energy  from  the  active  states  of  the 
surface . 

The  band  structure  of  an  Isolated  CO  overlayer  can  be 
found  in  a  previous  contribution  from  our  group  on  CO  bonding  to 
metal  surf  aces 1 2 ,  as  well  as  in  the  work  of  others  .  We  will 
focus  our  attention  on  the  coverage  and  site  dependence  of  CO 
chemisorbed  on  Rh(lOO).  The  results  of  our  calculations  are 

summarized  in  Table  l .  We  have  considered  two  adsorption 

geometries,  on -top  (3)  and  two -fold  bridging  (4),  although  it  is 
known  that  an  on  top  adsorption  is  strongly  preferred  for 
lower  coverages  on  the  Rh(lOO)  surface.  The  coverage  dependence 


analysis  was  carried  out  for  each  adsorption  geometry  with  three 
different  coverages  of  the  surface  of  l ,  1/2  and  1/4  which  is 
represented  by  a  top  view  in  5,  6  and  7,  respectively.  Table  1 

however  reports  the  results  for  one  coverage,  c(2x2)  and  0=0.5, 
because  the  results  for  other  coverages  are  so  similar. 
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Table  1. 


Some  Bonding  Characteristics  for  c(2x2)CO  on  a  Model 
Rli(  100)  Slab,  6:  0.5. 


CO  Electron 
Densities 


Overlap 

Populations 


Binding 
Energy  (eV)a 


a  Binding  energy 


Free  CO 

On  -top 

Bridging 

4c 

2 

1.75 

1 .72 

*"x 

2 

1.99 

1.95 

iXy 

2 

1.99 

1.95 

5 C 

2 

1.57 

1.59 

2S, 

0 

0.35 

0.45 

2Uy 

0 

0.35 

0.75 

C-Rli 

- 

0.57 

0.52 

C-O 

1.44 

1.23 

1.09 

2.45 

3.07 

defined  as  (Eco  ♦  ERh(lOO)>  ~  EC0+Rh(100)* 
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Q-  l.pflxl!  Q-0.5,c(2x2)  Q-025.p(2x2) 

on-top 

c(2x2)  p(2x1) 

two- fold 
bridging 

5  6  7 

Note  the  forward  donation  in  the  depopulation  of  the  5 a  (and  4c ) 
and  the  hack  donation  in  the  population  of  the  air*.  In  the 
bridging  mode  the  degeneracy  of  the  ir  orbitals  is ,  of  course , 
lost;  the  greater  back  donation  is  to  the  w  orbital  m  the 
Rh-C-Rh  plane,  which  is  able  to  interact  in  o  fashion  with  the 
Rh  orbitals.  With  that  greater  back  donation,  there  comes  a 
greater  weakening  of  the  C-0  bond.  The  binding  energies  cannot 
be  relied  on,  because  we  assume  and  cannot  optimize  a  Rh-C 

distance .  Not  much  changes  for  other  coverages ,  except  that  the 

binding  energy ,  by  its  decrease ,  clearly  shows  that  the  CO's  are 
uncomfortable  compressed  at  ©  =  l .  o . 
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Let's  look  with  a  little  more  detail  at  the  interactions 

in  the  on -top  c(2x2)CO+Rh( 100)  system.  At  left  m  Figure  1  is 

Figure  1  here 

the  DOS  (density  of  states)  of  a  naked  Rh(iOO)  slab  (Fig.  la),  m 
the  center  (Fig.  lb)  the  composite  surface  plus  CO  overlayer,  and 
at  right  (Fig.  lc)  the  energy  levels  of  an  isolated  CO  molecule. 

The  density  of  states  curve  clearly  shows  that  the  major  surface  - 
adsorbate  interactions  involve  5 o  and  2tt*  CO  orbitals,  although 
the  4o  orbital  is  also  shifted  toward  lower  energies.  The 
metal -CO  5 o  interaction  is  mainly  accomplished  through  electron 

donation  to  appropriate  surface  orbitals,  i.e.  d  which  are 
pointed  directly  toward  the  5 o  and  4o  orbitals  ( 0 ) .  The  net  result 
is  depopulation  of  the  5o  and  4c  from  the  2  electrons  each  carry 
for  a  free  CO  molecule ,  to  1.07  and  1 . 76  electrons ,  respectively , 
m  the  composite  system  (Table  1 ) .  As  a  result  of  this  the  5o  and 
4o  are  pushed  down  in  energy  (Figure  2a).  The  d^  orbitals 
are  pushed  up  and  consequently  a  greater  portion  of  the 


d  2  band  is  raised  above  the  Fermi  level.  Another  very  strong 

2 
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Energy  (eV) 


DOS - • 

a 


DOS  - - 

b 


c 


Figure  1.  (a)  Total  DOS  of  the  three  layer  slab,  with  Fermi 

level  (£f)  and  s,  p  and  d  bands  indicated;  (b)  Total  DOS  of  the 
c (2x2 ) CO-Rh ( 100)  system  (broken  line)  and  projected  DOS  of 
adsorbed  CO  molecule  (black  areas) ;  (c)  Energy  levels  of  an 
isolated  CO  molecule. 
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metal -adsorbate  interaction  is  through  the  co  an11  orbitals 

interaction  with  the  surface  dv  orbitals  (dxz,yz)>  as  represented 
schematically  in  9.  The  initially  empty  2ir”  orbitals  are  now 
partly  populated  in  the  composite  CO+Rh(iOO)  system,  with  an 
electron  density  of  0.76  electrons  .while  the  corresponding  surface 
d*  states  are  depopulated  . 

Pigures  2,3  here 

Another  way  to  follow  the  involvement  of  the  various 
orbitals  in  bond  formation  is  through  the  crystal  orbital  overlap 
population  (COOP)  analysis50.  The  COOP  curve  really  represents  an 
overlap  population  weighted  density  of  states.  It  weights  the 

states  in  each  energy  interval  by  their  contribution  to  the 
overlap  population.  Such  a  COOP  curve  for  C-o  and  C-Rh  bonding  in 
our  c  ( 2x2  )CO+Rh  ( 100)  system  is  presented  in  Figure  3.  The  main  CO 
orbital  contributions  are  marked  and  we  can  see  that  the  2w* 
orbitals  (i.e.  their  bonding  combination  with  the  dw)  are  c-Rh 
bonding  and  C-O  antibonding.  The  CO  is  orbitals  are  still  the 
major  bonding  CO  orbitals.  The  5 o  is  both  c-o  and  C-Rh  bonding, 
but  its  absolute  contribution  to  these  bonds  is  small.  On  the 

other  hand,  the  COOP  curve  shows  that  the  4o  makes  a  major 
contribution  to  C-Rh  bonding. 

What  is  the  reason  for  large  involvement  of  the  4o? 

One  thinks  4o  is  an  oxygen  based  lone  pair ,  5o  is  on 

carbon .  Actually  both  are  quite  delocalized ,  and  Table  2  shows 
their  composition  for  free  CO  with  our  parameters . 

Unfortunately,  the  make-up  of  these  orbitals  is  somewhat 
parameter  dependent .  One  can  see  the  delocalization  in  both . 
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Table  2 


C 


.  The  wave  function  coefficients  of  free  CO  4o  and  5 o 
molecular  orbitals  with  our  parameters  (see  Appendix) 


40 

50 

(S) 

-0.7197 

0.3306 

(px> 

0 . 0000 

0 . 0000 

(Py  ) 

0 . 0000 

0 . 0000 

(PZ) 

0.4051 

-0.4674 

(S) 

0.2949 

0.1147 

(PX) 

0 . 0000 

0 . 0000 

(Py) 

0 . 0000 

0 . 0000 

(PZ) 

0.4362 

0.7766 
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Figure  2.  Projected  DOS  (black  areas)  of  5a  (a)  and  2n  (b)  CO 
orbitals  in  a  c(2x2)CO-Rh(100)  system.  Broken  lines  represent 
the  total  DOS  of  the  system.  The  peak  between  -16  and  -17eV  is 
mainly  CO  4a. 
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Figure  3.  Crystal  Orbital  Overlap  Population  (COOP)  curve  for  an 
on-top  c ( 2x2 ) CO-Rh ( 100 )  system.  Full  line:  carbon-oxygen  bond; 
dotted  line:  rhodium-carbon  bonding. 
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Really  what  matters  is  the  overlap  these  orbitals  have 
with  another  group .  With  the  parameters  we  have  used ,  a  probe  H 
is  orbital  placed  at  a  distance  l .  0  -  2 .  o  A  away  from  the  carbon 
end  of  CO  actally  has  a  bigger  overlap  with  4a  than  5a .  Still  its 
net  interaction  with  4a ,  though  substantial ,  is  less  than  with  5o 
because  5o  is  closer  in  energy  to  the  hydrogen  (or  the  frontier 
orbitals  of  a  metal ) . 

The  CO  4a  orbital  contribution  to  chemisorption  in  the 
present  study  might  be  overestimated  by  the  use  of  charge 
iterated  parameters  for  carbon  and  oxygen  atoms  in  CO  overlayer 
on  the  rhodium  surface  (see  Appendix).  But  even  with  the  use  of 
standard  carbon  and  oxygen  parameters  in  extended  Huckel 
molecular  calculations51  we  had  almost  the  same  effect  as  here. 

We  note  here  also  that  others  have  suggested  a  significant  role 
for  the  4a  orbital 1 9a  •  * . 

CO  Bridging  Site 

We  consider  the  bridging  geometry  in  some  detail  not 
because  it  is  favored  for  CO  on  Rh  ( 1 00 ) ,  but  because  of  the 

possibility  that  such  a  state  may  be  Involved  In  the  coadsorption 
pattern  with  HO. 

Our  results  for  an  adsorption  geometry  presented  In  4  and 
for  a  c  ( 2x2 )  LEED  pattern  (lower  left  part  in  4)  are  also 
presented  in  Table  1 .  We  have  chosen  the  same  C  -Rh  distance 
although  it  is  Known  that  bridge -bonded  species  should  have 
somewhat  longer  distances . 

The  bonding  here  is  a  little  bit  more  complicated .  The  CO 
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2ir"  and  5 0  orbitals  still  interact  with  the  surface  d  2  orbitals 

z 

but  less  effectively .  The  corresponding  bonding  combinations  are 
shown  m  10  and  11.  In  addition ,  both  CO  frontier  orbitals  can 


interact  even  more  efficiently  with  the  surface  dXZ|yZ  orbitals; 
the  bonding  combinations  are  shown  in  12  and  13.  The  2**  orbitals 
interact  with  surface  antibonding  dXZ|yZ  orbitals  (13),  that  is 
with  the  upper  part  of  the  d  band.  Since  the  2*"  energy  falls  just 
at  that  energy,  the  interaction  is  quite  efficient.  Thus  the 

filling  of  the  co  2w"  orbitals  is  higher  in  the  bridging  that  in 
the  on -top  geometry  (Table  t).  Consequently,  the  C-0  overlap 
population  is  further  reduced  to  1.094  in  a  c(2x2)  pattern, 
following  the  well-established  trend  that  two -fold  bridge -bonded 
species  have  a  lower  stretching  frequency  than  terminally  bonded 
ones . 


12  13 
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It  is  very  interesting  to  note  that  the  4a  orbitals  are 

also  active  here.  Even  more  so  than  m  the  on -top  site;  note 
in  the  Table  1 ,  the  greater  depopulation  of  4a  for  CO  bridging . 

MO  on  the  RIM  100)  Surface 

The  chemisorption  of  nitric  oxide  on  transition  metal 
surf  aces  has  many  similarities  to  that  of  CO ,  but  is  more  complex 
due  to  the  additional  geometrical  possibilities  (linear , bent )  and 
the  additional  electron  occupying  the  2x“  orbital .  For  HO  on 


group  VIII 

transition 

metals  there  is  evidence 

for 

some  low 

temperature 

dissociative 

adsorption  and  a  variety 

of 

molecular 

bonding  geometries.  These  Include  a  form  with  the  NO  axis  tilted 
away  from  the  surface  normal  (bent ) ,  a  bridged  form  and  a  linear 


one . 

For 

Instance , 

the 

bent  HO  species 

were  found 

on 

reconstructed 

Pt(lOO)52 

and 

Hl(lll)53  surfaces. 

Both  bent 

and 

linear 

forms 

of  adsorbed  HO 

were  found  at  room 

temperature 

on 

polycrystalline  Ni30c. 

One  of  the  first  broad  investigations  of  HO  chemisorption 
was  carried  out  by  a .  Broden  et  al 1 .  According  to  the  line  drawn 
by  these  authors,  dividing  the  transition  metals  into  those  on 
which  adsorption  is  strictly  molecular  and  others  where 
dissociation  is  observed,  rhodium  (unmarked)  is  most  likely  to 
behave  as  Hi  or  Ir,  which  exhibit  both  molecular  and  dissociative 
HO  adsorption.  In  1979,  Castner  et  al.2  reported  that  HO 
chemisorption  on  the  Rh(lOO)  occurs  with  high  initial  sticking 
probability  at  25°C,  forming  a  c(2x2)  LEED  structure.  These 

workers  stated  that  thermal  HO  desorption  was  complex .  In  the 

same  year,  Campbell  and  White37  reported  that  a  polycrystalline 
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Rh  sample  at  saturation  has  about  15X  of  dissociated  HO . 

For  some  time  thereafter ,  MO  chemisorption ,  especially  on 
theRhdOO),  did  not  receive  much  attention.  Recently,  Ho 
and  White31  reported  a  low  temperature  (-100H)  adsorption 
and  decomposition  study  of  HO  on  Rh(iOO).  Adsorption  was 
predominately  molecular,  with  a  small  dissociative  contribution 
at  low  coverages,  and  a  saturation  coverage  of  0.65  ML.  The 
authors  suggested  that  there  are  two  major  NO  states  involved  in 
the  chemisorption.  The  first  corresponding  to  normal  molecular 
adsorption  and  the  second  has  highly  tilted  HO  molecules ,  such 
that  both  N  and  o  atoms  are  bound  to  the  surface  with  a  weakened 
H--o  bond.  Upon  heating  the  latter  decomposed  before  desorption. 
Later,  Villarrubla,  et  al.43  have  resolved  the  vibrational 
modes  of  920  and  1600  cm'1  for  two  such  different  HO  adsorption 
states  at  90K.  The  lower  frequency  state  was  attributed  to  a 
highly  inclined  or  side -on  (lying  down)  bonding  mode  of  HO.  A 
higher  frequency  state  of  about  1600  cm _1  (more  accurately  at 

1590  for  low  coverages,  shifting  to  1700  cm'1,  at  higher 

coverages )  was  not  unambiguously  assignable  to  a  specific 

adsorption  site  and  geometry .  Linearly  bonded  ho  stretching 
frequencies  in  transition  metal  nltrosyls  have  been  observed 
in  this  range,  although  most  reported  frequencies  on  single 
crystal  surfaces  are  slightly  higher .  On  the  other  hand ,  several 
previous  assignments  of  two -fold  bridged  adsorptions  are  close 
to  these  values ,  ranging  up  to  around  1 600cm  ~ 1 .  Very  recently , 
the  same  authors54  reported  a  broad  investigation  on  NO  adsorption, 
decomposition  and  desorption  on  the  RhdOO)  surface.  They  verified 
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the  existence  of  two  major  HO  states.  04  (lying  down)  and  a 2 
(vertically  bonded  one).  Upon  heating.  6 2.7.  of  HO  decompose, 

while  the  the  rest  desorb  molecularly  with  Ea=26i3  Kcal/mol.  The 
a |  HO  decompose  with  Ea  =  10. 510.7  Kcal/mol,  while  the  decomposition 
of  vertically  bonded  a  2  HO  is  more  complex,  strongly  coverage - 
dependant ,  and  goes  through  a  lying  down  a }  HO  state  as  an 
intermediate .  Por  a  more  detailed  discussion  of  this  subject  see , 
for  example,  references  32,  33,  43  and  54. 

It  appears  that  for  the  chemisorption  of  HO  on  the 
Rh(lOO)  surface  we  do  not  have  enough  experimental  information 
to  permit  an  unambiguous  assignment  of  adsorption  pattern  and 
geometry.  There  is  although,  some  evidence  that  HO  should  adsorb 
vertically  at  two -fold  bridging  positions,  since  Viliarrubia  and 
Ho54,  in  their  model  to  describe  HO  adsorption  on  Rh(iOO),  placed 
vertically  bonded  03  HO  in  the  two -fold  bridging  sites,  and  the 
lying  down  at  HO  in  the  four -fold  hollows.  Por  some  other  Rh 
surfaces,  such  as  Rh(lll),  Q.B.  Fisher  and  co workers 32  reported 
that  HO  chemisorbs  in  two -fold  bridge  sites  at  all  coverages  and 
that  HO  is  partly  dissociated  at  lower  coverages  (below  0.2 
monolayers)  with  an  activation  energy  for  dissociation  of  19 
Kcal/mol.  At  higher  coverages,  adsorption  is  molecular.  Therefore, 
in  order  to  examine  adsorption  and  bonding  of  linearly  bonded  HO  we 
performed  an  extended  Huckel  calculation  for  on -top  and  two -fold 
bridging  geometries,  for  several  coverages  ranging  from  9  :  l  to 
e  *  0.25.  The  adsorption  geometries  used  were  the  same  as 
for  CO  adsorption ,  as  indicated  in  3  and  4  and  with  three 
different  coverages  as  in  5,  0  and  7.  Table  3  chooses  again  one 
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Table  3.  some  Bonding  Characteristics  for  of  On -top  and 
Bridging  c(2x2)HO+Rh(  100)  system, 0:  0.5. 


Free  HO  On -top  Bridging 


NO  Electron 

40 

2 

1.55 

1 .61 

Densities 

iwx 

2 

1 .95 

1.96 

iSy 

2 

1.95 

1.59 

50 

2 

1.71 

1.73 

0.5 

1.09 

0.66 

2Sy 

0.5 

1.09 

1.27 

overlap 

H-Rfc 

- 

0.72 

0.54 

Populations 

N-O 

i  .  15 

0.69 

0.91 

Binding 
Energy  (eV)a 

- 

2.50 

3.24 

a  Binding  energy  deffined  as  (Bro  +  ®Rh(l00))  ~  ®HO+Rlx<  100)  • 
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typical  coverage  for  a  comparison . 

HO  chemisorption  on  Hi  ( 1 1  i )  has  been  examined  earlier  by 
this  group33.  The  molecular  orbitals  of  HO  are  very  much  URe 
those  of  CO  and ,  so  are  the  important  interactions .  Because  of 
the  extra  electron  and  the  low  energy  of  the  HO  as"  orbitals, 
the  total  filling  of  the  Zv*  derived  orbitals  of  the  composite 
system  is  much  higher  than  for  CO. 

This  system  is  less  coverage -dependent  than  the 
corresponding  CO  system .  This  is  indicated  by  calculations 
not  reported  here,  in  which  one  observes  smaller  changes  in  the  N-0 
overlap  populations  when  going  from  full  to  1/4  coverage. 
Experimentally ,  the  shift  m  stretching  frequencies  in  low -high 
coverage  systems  is  very  significant .  On  the  other  hand ,  it 
seems  that  our  prediction  as  to  which  of  these  two  different 
adsorption  geometries  is  preferred,  based  on  the  "binding 
energies"  .  is  in  agreement  with  the  experiment,  since  for 
Rh(iii)52  and  probably  Rh(tOO)5*,  the  HO  adsorption  goes 
exclusively  through  a  bridging  geometry .  For  this  geometry  our 
binding  energy  is  considerably  stronger . 

It  is  interesting  to  note  that  the  involvement  of  4 a  in 

adsorption ,  discussed  for  CO  chemisorption ,  is  also  present  here . 

The  total  effect  is  smaller  than  in  CO,  despite  the  smaller 

energy  separation  between  the  4 o  and  5a  levels.  As  for  CO, 

this  effect  is  stronger  in  the  case  of  bridging  geometry .  The 
4o  orbitals  are  more  depopulated  m  this  geometry  (Table  3). 
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The  Bending  of  MO 

Since  it  is  difficult ,  in  principle ,  to  resolve  experimentally  the 
problem  of  bridging  linear  versus  terminal  bent  NO  in  the  NO+Rh  (100) 
system ,  we  wanted  to  investigate  this  system  theoretically .  We  also 
wanted  to  study  the  proposed »  3 1  multi  -step  decomposition 
mechanism  of  vertically  bonded  NO ,  that  goes  through  the  NO  lying 
down  geometry  as  an  intermediate .  a  corresponding  analysis  of 
the  NO+Nl (ill)  system33,  and  of  the  bonding  in  discrete 
nitrosyls27  was  carried  out  by  this  group  earlier .  Thus  we 
emphasize  here  only  the  most  important  features . 

There  are  two  likely  ways  of  bending,  from  an  on -top  site 
and  from  the  bridging  position .  For  the  sake  of  simplicity ,  we 
did  all  the  calculations  for  the  (2x2)  geometry  so  that  the 
coverage  effects  while  bending  are  not  included. 


Bending  in  the  On -top  Site 

Let  us  consider  bending  in  an  on  -top  site  of  NO  adsorption , 
in  a  (2x2)  geometry  indicated  in  the  upper  drawing  m  7.  We  will 
assume  that  the  bending  takes  place  in  the  yz  plane  as  shown  in 
14  in  both  perspective  and  top  views .  The  bending  or  tilting 


14 
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angle  #  is  defined  between  the  surface  normal  and  the 


direction 


of  N-0  bond.  In  14  we  give  the  coordinate  system  that  will  stay 
unchanged  throughout  this  paper . 

The  results  of  our  calculations  are  shown  in  Table  4  as  a 

function  of  bending  angle  * .  As  the  bending  angle  increases  from 
0°  to  45°,  the  N-0  overlap  population  decreases  i.e.  the  NO  bond 
is  weakened .  The  N  -Rh  overlap  population  also  decreases .  The 
orientation  of  NO  5o  and  2tj  orbitals  and  the  dy2  and  d  ^  of  the 
surface  Rh  are  shown  in  15  and  16.  Upon  bending,  the  orientation 
of  the  5o  orbital  enables  it  to  interact  with  metal  yz  orbitals , 

but  this  interaction  is  not  strong  enough  to  compensate  for  the 
loss  in  5o  -d  3  interaction  ( 6 ) . 


a 


15 


b 


16 


As  the  NO  bends,  4 o  and  5o  interact  less  with  the  metal. 

This  is  why  their  population  moves  toward  2.0.  The  wy  and  o 

orbitals  mix  and  hybridize,  so  that  2iry  (really  a  mixture  of  2ny 

and  2vz )  interacts  with  more  orbitals  of  the  metal ,  such  as 

d  2 »  a  and  pz .  Therefore ,  it  gains  more  electron  density  than 
z 

the  2«x  orbital .  This  increased  occupation  of  antibonding 
levels  is  the  main  reason  for  weakening  of  the  N-0  bond  in  the 
bent  geometry . 
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Table  4.  Bending  of  HO  in  the  yz  plane  (along  the  surface  Rh-Rh 
bond)  in  an  on -top  site  (see  14). 


* 

0o 

45° 

90° 

HO  Electron 

49 

1.95 

1.89 

1 .90 

Densities 

mx 

1.99 

1.98 

1.97 

lwy  (z) 

1.96 

1.98 

1.90 

59 

1.71 

1.77 

1 .83 

2*X 

1.09 

1.02 

1.37 

2wy  (z ) 

1.09 

1.43 

1.09 

overlap 

H-O 

0.99 

0.82 

0.85 

Populations 

H-Rh(l) 

0.72 

0.94 

0.42 

0  -Rh  ( 2 ) 

- 

0.00 

0.12 

dO-Rh(2) 

4.07 

3.30 

2.44 

Binding 

2.50 

2.21 

1.52 

Energy  (eV ) 
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17 

By  further  bending ,  +  -  90° ,  we  reach  the  situation  shown  in 
17.  The  NO  bond  now  lies  parallel  to  the  Rh(l)-Rh(2)  surface 

bond ,  with  the  N  atom  being  exactly  on  -top  of  the  Rh  ( l )  atom  . 
This  geometry  is  particularly  important  since  some  similar  states 
have  been  detected  experimentally .  This  will  be  discussed  in 
more  detail  below ,  in  a  section  analyzing  possible  geometries  of 
NO  parallel  to  the  surface .  Upon  this  further  bending  ,  the  N  -O 
overlap  population  increases  somewhat,  the  N-Rh  overlap 
population  decreases  further.  On  the  other  hand,  the  0-Rh(2) 
overlap  population  grows  to  the  substantial  positive  value  of 
0.12  with  the  relevant  atoms  being  2 . 44  A  apart .  Now  we  have  a 
very  interesting  bonding  pattern  in  which  the  2wz  interacts 
with  the  surface  in  the  way  the  5 a  did  in  the  linear  geometry , 
and  the  5 a  takes  over  the  interaction  with  the  Rh  ( l )  dyz 

orbitals.  The  total  NO -surface  interaction  is  weaker  and  the 
binding  energy  is  reduced  to  only  1.52  (Table  4).  We  also 
studied  further  bending  of  NO  up  to  110°.  The  results  are  not 
reported  in  the  table  but  they  continue  the  trend  of  N-0  and  N-Rh 
overlap  populations  decreasing  further.  It  appears  that  this 

bending  is  not  a  good  candidate  for  a  dissociation  mechanism 
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pathway  since  the  changes  in  the  N-0  overlap  population  do  not 
show  the  required  significant  NO  bond  weakening . 

Another  possibility  for  bending  In  an  on -top  position  is 
shown  in  16.  Now  the  bending  is  toward  a  four -fold  hollow,  l.e. 
m  the  diagonal  of  the  xz  and  yz  planes .  One  is  led  to  consider 

this  variant ,  as  a  precursor  to  N  -O  dissociation  ,  because  it  is 
known  that  atomic  species  tend  to  occupy  high  symmetry  positions 
such  as  the  four-fold  hollow  site55. 


top  view 

18 

What  happens  in  the  calculations  is  that  the  binding  energy 
decreases  steadily  and  the  NO  overlap  population  behaves  Just  as 
it  did  for  bending  "along  a  bond”.  So  this  bending  toward  a  hollow 
does  not  look  like  a  good  prospect  for  the  dissociative  reaction 

pathway  either . 

Bending  in  a  Bridging  Position 

The  two  bending  modes  we  single  out  for  study  are  in  the  yz 
plane  (toward  a  hollow,  19)  and  in  the  xz  plane  (toward  a 
bond,  20).  The  first  process,  bending  toward  a  hollow,  does  not 
seem  to  lead  to  anything  new  -  once  again  the  NO  bond  strengthens 
on  90°  bending. 
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The  second  bending  mode ,  20 ,  is  more  interesting; 


the 


19  20 

computational  results  for  two  bending  angles  are  summarized  in 
Table  5.  There  are  two  entries  for  the  N-Rh  overlap  populations 
since  the  bending  is  not  symmetrical  with  respect  to  the  the 
Rh  ( i  )  and  Rh  ( 2 )  surface  atoms .  As  NO  bends ,  the  orientation  of 
the  o  and  s  orbitals  changes .  Among  these ,  the  2vy  stays 
essentially  the  same,  but  there  is  a  certain  drop  in  the  2irx 
(i.e,  2*x,z  for  0 s45® )  electron  density.  As  a  result,  the  N-0 

overlap  population  becomes  larger.  The  bending  of  NO  in  the  x 

direction  reduces  the  possibility  of  H-Rh(2)  interaction  since  5 o 
is  moved  away  and  directed  more  toward  Rh  ( 1 ) .  The  corresponding 
overlap  population  is  significantly  reduced.  The  N-Rh(l)  overlap 
population  is  somewhat  increased.  The  overlap  population  between 
oxygen  and  the  Rh(2)  atom  is  still  negative  so  that  we  still  have 
an  antibonding  interaction  here . 

With  further  bending  of  NO  to  90°,  the  situation  changes  a 
lot.  The  4o-5o  interactions  with  the  surface  Rh(i)  and  Rh(2) 
orbitals  are  increased  and  both  of  these  orbitals  are 
depopulated .  There  is  a  certain  increase  in  lir  orbital 
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Table  5 


Bending  of  NO  m  the  xz  plane  while  bridging  (20). 


0° 

45° 

90° 

NO  Electron 

40 

1.51 

1.83 

1 .73 

Densities 

l*X  (Z) 

1.89 

1.67 

1.88 

lXy 

1.95 

1.95 

1.81 

50 

1.73 

1.75 

1.64 

2wX(Z) 

1.26 

0.95 

1.27 

2l?y 

0.67 

0.66 

0.94 

Overlap 

N-0 

0.91 

1.01 

0.89 

Populations 

N-Rh(l) 

0.54 

0.56 

0.49 

N  -Rh  ( 2  ) 

0.54 

0.32 

-0.05 

O-Rh  { 2 ) 

- 

-0.03 

0.71 

<*0-Rh(2) 

2.83 

2.22 

1 .37 

Binding 

3.21 

2.67 

3.05 

Energy  (ev ) 
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interactions .  The  2iry  orbitals  are  still  almost  unchanged ,  but 
the  2*x  (the  2sz  now)  occupancy  is  getting  larger.  The  2sz 
interacts  in  the  way  the  5a  interacted  in  the  linear  case,  and  is 
filled  a^rell .  The  net  result  of  all  of  this  is  a  significant 
weakening  of  the  NO  bond  (note  the  decrease  in  the  NO  overlap 
population).  N-Rh(l)  overlap  population  decreases.  There  is  a  very 
weak  antibonding  N  -Rh  ( 2 )  interaction ,  but  a  large  increase  in  the 
O-Rh ( 2 )  overlap  population,  occurring  mainly  through  the  4a  and 
surface  dxzmteractlons .  It  must  be  said  that  the  O-Rh (2)  distance 
is  much  too  short  at  90° ,  so  that  bonding  is  overestimated . 

The  binding  energy  increases  at  90°  in  mode  20,  and  finally 
we  have  found  a  way  to  weaken  the  NO  bond .  Perhaps  this  is  the 
way  the  NO  dissociates .  But  first  let's  examine  a  set  of 
alternative  geometries . 

LTlPM  Down  BO 

In  the  preceding  discussion ,  we  have  mentioned  a  few 
different  geometries  with  the  NO  lying  down,  l.e.  parallel  to 
the  surface.  Such  NO  species  stable  at  lower  temperatures  and 
coverages,  have  been  observed  experimentally  on  this  surface. 
First,  Ho  and  White31  proposed  an  NO  adsorption  state  with  a  weak 
NO  bond  and  with  both  N  and  O  interacting  with  the  surface  atoms 
for  several  desorption  peaks  observed  in  the  thermal  desorption 
of  NO  from  Rh(lOO).  Villarrubla,  et  al.**«5*  recently  reported 
a  stable  low  frequency  NO  species  adsorbed  on  Rh(lOO)  and 
assigned  it  to  either  a  highly  inclined  structure  or  one  lying 
parallel  to  the  surface,  and  proposed  a  model  where  lying  down 

NOs  are  occupying  four -fold  hollow  positions.  A  similar  Ng 
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adsorbed  species  have  been  observed  on  Fe  ( i  i  l ) 55  and  on 

X.  -promoted  Ru(OOl)57.  Side -on  bonded  CO  molecules  have  also  been 
reported  on  the  Cr(llO)  surface5®.  Such  experimental  findings 
have  been  theoretically  considered  in  the  case  of  CO  on  the 
Cr(UO)  surface  by  Hehandry  and  Anderson15  and  in  the  case  of 
HO  on  various  platinum  surfaces  by  the  Hasel  group59.  The  first 
have  used  an  atom  superposition  and  electron  delocalization  HO 
theory ,  with  the  surface  represented  by  a  two -layer  thick 
cluster  of  33  chromium  atoms. 

If  we  begin  to  study  the  bonding  possibilities  for  HO  lying 

down  on  a  surface,  we  are  led  naturally  to  consider  a  set  of 
geometries,  21-20.  21,  22  and  24  have  been  discussed  in  the 

previous  section .  These  include  not  only  fpolnt  reached  at  90° 
from  on -top  or  bridging  linear  origins,  but  other  geometries  as 
well . 

Some  computational  results  for  these  geometries  are  compared 
in  Table  6 .  Let's  discuss  them  one  by  one.  In  21,  obtained  by 
bending  from  the  linear  geometry,  we  lose  the  H -surface 
interaction  (the  H-Rh  overlap  population  gets  smaller)  and  gaining 
some  new  oxygen -surface  interaction  (see  Table  4) .  As  we  bend  HO 
from  linear  to  lying  down  position  22,  starting  in  the 
bridging  position,  the  nitrogen  completely  loses  its  bonding 
interaction  with  the  rhodium  atom  toward  which  the  HO  bends, 
the  oxygen  atom  comes  very  close  to  this  Rh  atom ,  acquiring  a 
large  positive  overlap  population  (see  Table  5).  It  seems  likely 
that  the  HO  group  as  a  whole  will  shift,  in  response  to  bonding 
effects,  to  a  geometry  such  as  23.  For  this  geometry,  we  have 
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Table  6.  various  possibilities  of  NO  lying  parallel  to  the 


Rh(l00)  surface. 


Geometry 

overlap 

Populations 

do-Rn(A> 

Binding 
Energy  (eV ) 

N-0 

N-Rn 

O-Rh 

21 

0.85 

0.42 

0.12 

2.44 

1.52 

22 

0.89 

0.49 

0.71 

1.37 

3.05 

23 

0.90 

0.55 

0.39 

1.90 

2.49 

24 

0.85 

0.40 

0.02 

2.79 

1.39 

25 

0.98 

0.52 

0.38 

1.90 

2.09 

25 

0.96 

0.34 

0.07 

2.44 

0.97 

27 

0.94 

0.47 

0.31 

1.90 

1.93 

25 

1.00 

0.24 

0. 12 

1.90 

1.98 
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Another  bending  possibility  from  an  on -top  position  is 
toward  a  four -fold  hollow,  leading  to  24.  The  bending  is  also 
accompanied  by  a  decrease  in  N-Rh  and  N-0  overlap  populations. 
This  bending  is  not  energetically  favored ,  at  least  not  if  the 
position  of  the  H  atom  is  fixed  during  the  bending .  one  possible 
change  might  be  to  move  the  H  atom  along  the  diagonal  of  the 
(100)  face,  to  reach  the  geometry  shown  in  25.  In  this  case  too, 
the  H-Rh  and  O-Rh  distances  are  i.90  A,  and  the  N-0  distance  is 
1.15  A  as  before . 

It  is  interesting  to  note  here  that  according  to  the 
calculations  of  Kehandru.  and  Anderson 1 5 ,  the  most  stable  form  of 
CO  adsorbed  on  Cr(UO)  is  with  CO  lying  across  the  longer  diagonal 
of  the  (110)  face.  Similarly  the  Masel  group59  has  found  that 
there  is  a  good  match  between  the  orbitals  of  the  metal  and  of 
the  adsorbed  NO  (in  the  case  of  Pt(lOO))  if  the  NO  lies  across 
the  diagonal  of  the  unit  cell,  and  concluded  that  this  will  lead 
to  NO  dissociation.  In  going  from  24  to  25,  the  2w* -surface 
interaction  decreases  but  all  the  other  NO  orbitals  are  more 

involved  in  bonding ,  especially  the  4c  and  5c .  As  a 
consequence  of  this,  the  N-0  overlap  population  is  much  larger 
in  25  and  so  are  the  N-Rh  and  O-Rh  overlap  populations. 

This  movement  is  energetically  favored . 

Another  possible  geometry  considered  for  lying  down  NO 
originates  from  bending  from  the  bridging  position  toward  the 
four -fold  hollow  19.  The  top  view  of  HO  lying  down  in  this 
geometry  is  shown  in  25 .  During  this  bending ,  the  binding  energy 
increases  continuously ,  the  N  -O  overlap  population  increases 
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after  an  initial  drop  at  j*:45°  and  the  H-Rh  overlap  population 
constantly  decreases.  Since  by  this  bending  we  are  gaming  some 
new  attractive  O -surface  mteactions,  one  again  expects  further 
motion  to  optimize  these  interactions .  The  HO  group  can  move 
"right”  from  the  position  shown  m  26  to  reach  a  more  symmetrical 
position  .that  shown  in  27,  or  it  can  move  "left"  to  occupy  a 
position  where  it  lies  perpendicular  to  the  Rh  -Rh  surface  bond , 
directly  above  its  midpoint.  as  shown  in  26.  in  going  from 
26  to  27,  the  interaction  of  4o  and  5o  with  the  surface  orbitals 
increases ,  the  is  orbital's  electron  densities  decrease ,  while  the 
2ii”  increase  slightly.  As  a  result  of  these  changes,  the  N-0 

overlap  population  decreases  somewhat  and  overlap  populations 
between  oxygen  and  nitrogen  atoms  with  the  surface  increase.  The 
H-Rh  and  O-Rh  overlap  populations  in  Table  6  refer  to  the  nearest 
surface  Rh  atoms.  The  changes  in  these  overlap  populations  lead  to 
a  higher  binding  energy  and  this  movement  is  energetically  favored . 
The  same  is  true  on  going  from  26  to  26.  But  the  bonding  pattern 
and  consequently  the  ho  o  orbitals  densities  are  quite  different . 

The  HO  o  orbitals  do  not  have  a  suitable  geometry  for  interaction 
with  the  surface  atom  any  more.  The  overall  interaction  goes 
primarily  through  HO  «  orbitals,  but  the  total  HO  2v*  orbital 
filling  is  slightly  smaller  that  in  26  and  27,  so  that  the  N-0 

overlap  population  is  higher.  still  the  total  energy  change 
is  essentially  the  same  as  in  27. 

In  conclusion  we  can  say  that  a  geometry  for  highly  inclined 
or  lying  down  HO  on  the  Rh(iOO)  surface,  according  to  our 
calculations,  is  likely  to  be  23,  although  25  and  27  should 
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not  be  disregarded  as  possible  candidates .  It  should  be  noted 
again  that  27  was  proposed54  to  be  the  lying  down  HO  geometry.lt 
is  very  hard  to  give  a  definitive  answer  to  this  question  using 
our  computational  method .  The  discussion  in  this  section  was 
partly  based  on  energetic  arguments .  In  general  we  tend  to 
avoid  such  because  of  our  unwillingness  to  trust  the  total 
energies  given  by  the  extended  Huckel  method . 

The  Coadsorption  of  IQ  and  CO 

The  oxidation  of  CO  and  simultaneous  reduction  of  NO  over 
platinum  group  metal  surfaces  is  of  great  importance  in 
pollution  control  since  NO  is  a  major  pollutant .  There  have  been 
relatively  few  mechanistic  studies  examining  the  NO/CO  process 
on  well -characterized  rhodium  surfaces,  although  this  metal  is 
the  only  catalyst  which  works  effectively  to  reduce  NO  m 
automobile  exhaust  environment. 

Campbell  and  White37  have  examined  nitric  oxide  reduction 
with  carbon  monoxide  on  polycrystalline  rhodium  wire,  using  the 
flash  desorption  spectroscopy  -FDS .  They  found  that  adsorbed  NO 

molecules  dissociate  into  adsorbed  nitrogen  and  oxygen.  Oxygen 

reacts  with  either  adsorbed  or  gas -phase  CO  to  form  carbon 
dioxide.  Dubois  et  al.5  have  reported  evidence  for  an  oxygen 
intermediate  in  the  NO/CO  process  on  the  Rh(331)  single  crystal 
surface  using  HRBBLS  and  AES  and  indicated  that  under  their 
reaction  conditions,  the  formation  of  gaseous  molecular  nitrogen 
is  a  fast  process .  Two  recent  mechanistic  studies  came 
from  the  work  of  3.B.  Fisher 35b  •  35c  » 35e  and  Hendershot  and 
Hansen 3® .  The  first  studied  the  reaction  on  the  Rh  ( i  1  i ) 
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surface  and  concluded  that  the  rate  determining  step  (as  wen 
for  co/Og  reaction)  is  co (a)  ♦  O(a)  -->  co2(gas).  in  recent 

worK36e  they  argued  about  the  possibility  of  true  bimolecuiar 
reaction  being  involved  on  Rh  ( l  i  l )  surface ,  and  that  the 
reaction  rate  on  Rh(lO0)  is  significantly  different  from 
the  one  on  Rh  ( 1 1 1 )  surface .  The  second  contribution  is  more 
interesting  for  us  since  it  deals  with  CO/MO  processes  on  the 
Rh(lOO)  surface.  The  authors  found  that  this  process  is 
selective  toward  the  production  of  Ng  and  a  COg,  but  that  some 
other  steps  must  be  present  in  the  mechanism ,  including  the 
formation  and  decomposition  of  NgO  and  HOg  adsorbed  species.  Also, 
they  stated  that  the  NO  dissociation ,  i.e.  adsorbed  atomic  oxygen 
does  not  play  an  important  role  in  the  high  pressure  reactant 
partial  pressure  conditions.  Vlllarrubla  and  Ho54  also  reported 
a  low  coverage  coadsorption  effects  of  lying  down  NO  with  Og ,  NO 
and  CO,  and  concluded  that  coadsorbates  will  convert  the  lying 
down  aj  NO  into  a  vertically  bonded  one,  with  CO  being  about  two 
times  less  effective  than  Og  and  NO. 

Now  let  us  investigate  the  coadsorption  of  CO  and  NO  on  the 
Rh(ioo)  surface  with  our  calculational  tools.  The  general 
procedure  that  we  will  adopt  is  to  put  CO  and  NO  on  the  surface 
at  the  same  time  in  representative  positions  suggested  by 
experimental  studies .  By  following  the  trends  in  orbital 
electron  density  and  overlap  population  changes,  we  can  look  for 
indicators  of  NO  bond  weakening  and  COg  formation . 

There  are  a  large  number  of  possible  arrangements  of  CO  and 
NO  on  this  surface ,  as  both  can ,  in  principle ,  assume  bridging 
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or  on -top  adsorption  sites  .  Let  us  first  consider  the  situation 
where  CO  and  NO  are  3.80  A  apart  (the  surface  unit  cell  diagonal). 
Two  possibilities  are  one  with  both  CO  and  NO  being  m  bridging 
positions  29  and  one  with  both  in  on -top  positions  30.  The 

over  layer  has  a  c(2x2)  structure  in  both  cases  comprised  of 

P  ( 2x  2 )  structures  of  the  individual  adsorbates  giving  a  total 
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coverage  of  1/2.  The  drawings  give  a  perspective  view  showing 
only  the  surface  layer  of  the  Rh  slab.  The  big  circles  represent 
Rh  atoms ,  while  c ,  N  and  O  atoms  are  indicated .  The  coordinate 
system  is  the  same  that  was  shown  in  14. 

In  these  geometries,  CO  and  NO  do  not  "feel”  each  other  very 

much .  The  overlap  populations  and  electron  densities  have  not 
changed  significantly  from  the  case  of  individual  adsorption. 
This  is  consistent  with  our  previous  conclusion  that  the 
coadsorbate  interactions  for  both  the  CO  and  NO  are  not 
significant  for  distances  greater  than  the  van  der  waais  radii . 

For  29  we  considered  the  bending  of  NO  with  the  bending 
angle  0  being  defined  as  before.  The  bending  direction  is  toward 
the  four -fold  hollow,  so  that  Cqo  and  Oho  would  come  closer  as 
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the  bending  angle  approaches  90°. 

With  bending .  the  NO  and  N  -Rh  overlap  populations  do  not 
seem  to  be  affected  much  by  the  presence  of  CO ,  since  they  are 
nearly  the  same  as  they  were  in  19.  The  C-0  overlap  population 
increases  and  the  C-Rh  overlap  population  decreases  relative  to 
the  case  of  CO  adsorption .  There  is  a  very  weak  repulsive 
interaction  between  the  adsorbates .  The  differences  are  small 
compared  to  1/2  coverage  adsorption  of  each  individual  adsorbate . 

The  origin  of  this  so-called  "repulsion”  may  be  competition  of 
the  CO  and  NO  for  surface  electron  density . 

A  similar  situation  is  observed  when  both  CO  and  NO  are  in 
on -top  positions  30  with  the  bending  direction  of  NO  being  along 


30 

the  surface  Rh-Rh  bond. 

It  is  clear  that  in  these  geometries  with  a  total  coverage 

of  1/2 ,  CO  and  NO  are  too  far  apart  to  interact  strongly . 
Therefore ,  to  model  the  reactivity  of  CO  and  NO ,  these  two 
adsorbates  must  approach  each  other  in  a  realistic  way.  One 
possibility  has  the  same  initial  geometry  as  in  30  but  this  time 
the  bending  is  toward  the  four  -fold  hollow  and  ,  hence ,  the 
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oxygen  of  the  NO  moves  closer  to  the  CO  molecule .  For  a 
bending  angle  of  90°,  the  Ono~cco  distance  becomes  2.65A, 
and  the  overlap  population  between  them  is  still  negligible . 

We  can  see  some  of  the  coadsorption  effects ,  but  they  are  still 
weak.  The  adsorbates  may  be  too  far  away  for  a  significant  through 
-space  interaction  or  this  geometry  may  not  be  suitable  for 
effective  orbital  interaction.  To  achieve  the  values  of  O-C-O 

angle  in  the  transition  metal  complexes  of  CO  2.  the  NO  bending 

angle  was  increased  to  110°.  No  significant  changes  were 
produced ,  however . 

Several  other  possibilities  exist  which  will  provide  for 
shorter  contacts  between  the  Cco  and  Ogo»  These  include 
increasing  the  total  adsorbate  coverage  or  making  adjustments  in 
the  adlayer  structure.  The  previous  discussion  was  based  on  an 
overall  c(2x2)  pattern  comprised  of  p(2x2)  patterns  of  each 
adsorbate . 
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We  can  imagine  such  a  situation  where  CO  is  adsorbed  on -top 
and  NO  m  the  two -fold  bridging  position  as  shown  in  31.  These 

are  the  sites  that  CO  and  NO  would  prefer  in  the  absence  of 
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competing  coadsorbate  interactions .  The  results  of  our 
calculations,  for  this  geometry  are  summarized  in  Table  7.  The 

coadsorption  effect  for  #  s  0°  is  somewhat  stronger  now,  and  is 
reflected  m  higher  c-0  and  N-0  overlap  populations.  This  is 
consistent  with  HREELS  data  showing  higher  CO  and  KO  stretching 
frequencies  with  increasing  coverages,  i.e.,  with  the  presence  of 
repulsive  interactions.  Some  decrease  in  C -surface  and  H -surface 
interaction  is  expected,  but  in  our  calculations,  these  overlap 
populations  are  essentially  the  same.  The  H-Rh(3)  overlap 
population  is  somewhat  higher  than  when  HO  is  adsorbed  alone  and 
is  suspected  to  be  a  consequence  of  CO -MO  repulsion.  With 
further  bending,  we  gain  some  attractive  C-Ogo  interaction  but 
the  repulsive  interactions  are  stronger,  so  that  the  C -surface 
interaction  is  reduced  and  C-0  overlap  population  increased.  The 
energetics  of  these  changes  show  a  constant  decrease  in  binding 
energies,  so  that  we  could  conclude  that  this  arrangement 
does  not  lead  to  desired  changes ,  either .  These  coadsorption 
effects  can  support  the  experimental  finding  of  Viliarrubia  and 
Ho54,  that  lying  down  HO  will  be  converted  into  a  vertically 
bonded  one  when  coadsorbed  with  CO. 

Another  possible  arrangement  of  coadsorbed  CO  and  HO  is 
shown  in  32.  Both  adsorbates  are  in  bridging  positions,  but  now 
the  overall  adsorption  geometry  gives  a  2x1  structure  with 
nearest  neighbor  coadsorbate  interactions  of  2.69  A.  The  HO 
bending  direction  is  toward  the  four  -fold  hollow . 

Again  when  HO  is  linear  there  is  no  significant  coadsorption 
effect.  CO  shows  the  same  effects  whether  the  nearest  neighbor 
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Table  7.  The  coadsorption  of  CO  and  HO  on  the  Rh(ioo)  surface 


<3i) 

with  the  HO 

bending 

toward  the 

four  -fold 

hollow . 

* 

oO 

45° 

'  90° 

110° 

CO  Electron 

50 

1 .67 

1.67 

1.64 

1.65 

Densities 

0.35 

0.36 

0.32 

0.66a 

ZVy 

0.36 

0.37 

0.34 

HO  Electron 

50 

1.72 

1.75 

1.75 

1.79 

Densities 

2lTX(Z) 

0.65 

0.99 

1.04 

1 .98a 

2l!y 

1.27 

1.29 

1.06 

Overlap 

c-o 

1.24 

1.23 

1.25 

1.26 

Populations 

C-Rh(l) 

0.66 

0.85 

0.61 

0.63 

H-O 

0.92 

0.69 

0.96 

1.00 

H-Rh  ( 3 ) 

0.55 

0.50 

0.34 

0.27 

C-°ho 

0.00 

0.03 

0.05 

0.01 

Binding 

6.26 

5.60 

3.43 

3.26 

Energy  (eV) 


a  The  total  over  two  Sir  orbitals . 
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adsorbate  is  NO  or  another  CO .  Of  course ,  things  change  when  NO 
starts  to  bend  since  the  C  -Oho  distance  is  getting  smaller . 


32 

For  the  0  =0°  the  distance  between  carbon  and  nitric  oxide  oxygen 
is  2.92  A  and  when  NO  bends  to  45° ,  it  is  reduced  to  2.02  A .  At 
this  distance ,  we  can  expect  some  interaction ,  and  it  is 
reflected  in  the  small  positive  C-Ojjo  overlap  population  of 
0.02.  The  C-o  and  C-Rh  overlap  populations  are  reduced.  Also, 
there  is  a  small  repulsive  interaction  between  the  carbon  and 
nitrogen  (a  small  negative  overlap  population  of  -0.02).  When 
NO  is  lying  down ,  the  carbon  and  nitric  oxide  oxygen  are  now 
only  1.54  A  apart  and  the  orientation  of  the  CO  irx  and  NO  5o 
orbitals  is  suitable  for  a  strong  interaction .  This  results 
in  a  large  positive  C-Oro  overlaP  population  of  0.45,  compared 
to  the  value  for  a  C-0  single  bond  of  0.55.  The  occupation  of 
the  CO  2«*  orbitals  increases  and  the  bond  strength  decreases . 

The  C  -surface  interaction  is  weaXer,  since  the  CO  2*”  orbitals 
now  interact  with  NO .  But ,  as  in  the  case  when  NO  is  lying  down 
alone  on  the  surface,  the  N -surface  interaction  is  highly  reduced 
and  even  additional  interaction  with  CO  is  not  strong  enough  to 
weaken  the  NO  bond  further.  A  significant  bonding  interaction  is 


Page  44 


observed  between  Cgo  and  0N0  •  though . 

Another  possibility  of  CO  and  HO  coadsorption  Is  shown  in 
33.  Here  both  CO  and  HO  are  in  on -top  positions  on  two  nearest 
surface  Rh  atoms.  The  CO -HO  distances  are  the  same  as  in 
previously  discussed  geometry  32 .  Again .  the  HO  bending  is  in  the 
xz  plane  and  the  calculations  were  performed  for  the  three 
bending  angles .  The  results  are  summarized  in  Table  6 . 


This  situation  is  reasonable,  as  CO  prefers  on -top  geometries 
on  Rh(100),  and  on  the  Rh(lll)  surface  HO  shifts  from  two -fold 
bridging  positions  to  on -top  sites  when  oxygen  is  present  on  the 
surface39. 

Let  us  first  discuss  the  coadsorption  effects  for  #:0°  and 

compare  the  results  to  those  for  individual  adsorption  of  CO  and 

HO.  There  are  no  changes  in  4c  and  li  orbitals  of  CO  and  HO 
(not  shown  in  the  Table).  In  this  geometry,  the  degeneracy  of  the 
2x*  orbitals  is  removed  and  these  are  considered  separately.  The 
CO  2e*  are  less  populated  than  for  CO  on  Rh(iOO)  (see  Table  l) 
and,  consequently,  the  C-0  overlap  population  is  larger.  The 
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Table  ©.The  coadsorption  of  CO  and  KO  on  the  Rh(ioo)  surface  in 


the  on -top  positions  with  NO  bending  in  tbe  xz  plane  (see  33) . 


0<5 

45° 

90° 

110° 

CO  Electron 

50 

1.66 

1 .67 

1.64 

1.61 

Densities 

2»X 

0.33 

0.43 

0.73 

0.59 

2Wy 

0.35 

0.35 

0.42 

0.36 

NO  Electron 

50 

1.70 

1.72 

1.55 

1.62 

Densities 

2irx(z) 

0.99 

1.23 

0.96 

1.02 

2*y 

0.97 

0.96 

1.25 

1  .48 

Overlap 

C-0 

1.24 

1.20 

1.02 

1  .  12 

Populations 

C-Rh(l) 

0.66 

0.67 

0.69 

0.6© 

H-0 

0.91 

0.67 

0.66 

0.82 

N-Rh  ( 2 ) 

0.72 

0.66 

0.42 

0.26 

C-N 

0.01 

-0.03 

-0.05 

-0.04 

c-oHo 

0.00 

0.04 

0.47 

0.38 

oHO-Rh ( i ) 

0.00 

0.05 

-0.03 

0.00 

Binding 

5.34 

3.61 

2.50 

2.96 

Energy  (eV) 
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same  is  true  for  HO  (see  Table  3).  Both  the  C-Rh  and  H-Rh 
overlap  populations  are  smaller .  This  represents  a  weaker 
interaction  with  the  surface  and  would  result  in  a  greater  CO  or 
HO  stretching  frequency. 

The  coadsorbate  interaction  is  stronger  when  HO  bends  to 
45°.  The  C-0  overlap  population  is  significantly  reduced  due  to 
the  Cco~Oho  interaction  which  now  has  a  positive  overlap 
population  of  0.04.  The  H-Rh  and  C-Rh  overlap  populations  are 
reduced ,  but  the  H  -o  overlap  population  is  greater  than  in  the 
case  of  single  HO  adsorption .  Evidently ,  these  changes  are  a 
consequence  of  overall  repulsive  interactions . 

With  further  bending  of  HO  to  90°,  the  C-0  and  C-Rh  overlap 
populations  are  reduced  further.  The  H-0  overlap  population  is 

somewhat  larger  and,  in  fact,  greater  than  when  the  HO  is 

adsorbed  alone  in  the  same  geometry  (0.60  vs.  0.65).  We  have 
significant  interaction  between  carbon  and  nitric  oxide  oxygen 
with  an  overlap  population  of  0.47.  The  Cco_0HO  distance  is  the 
same  as  in  32 ,  but  here  they  interact  more . 

For  #s0°-90°,  the  trend  m  binding  energy  is  toward  weaker 
adsorption,  but  with  further  bending  to  liO°,  the  overall 
interaction  becomes  more  attractive ,  with  binding  energy  rising 
again.  The  C-o  overlap  population  is  somewhat  larger  but  the 
C-Rh  interaction  is  reduced  further.  The  H-0  overlap  population 
is  smaller  (0.62)  than  when  the  HO  is  alone  at  the  surface  m 
the  same  position  (0.64). 

In  general ,  with  this  bending ,  we  are  approaching  the 
situation  where  there  is  a  tendency  for  CO 2  formation,  but  the 
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energetic  changes  In  these  systems  are  not  favorable .  We  even 

tried  to  calculate  this  system  with  smaller  initial  bending 
angles  of  1 5°  and  30°  to  investigate  the  energy  changes ,  but  the 

energy  difference  of  these  systems  (5.23  and  4.65  ev , 
respectively)  indicate  a  constant  decrease  in  binding  energies 
with  bending . 

We  can  also  imagine  coadsorption  geometries  where  both 
the  CO  and  HO  are  in  the  bridging  positions,  being  bound  to  the 
same  surface  Rh  atom .  one  of  two  such  possibilities  is  shown  in 
34,  where  CO  and  HO  are  bridging  mutually  perpendicular  surface 
Rh-Rh  bonds.  Here  we  can  imagine  two  possible  ways  of  HO  bending 


34 

in  the  xz  plane  and  the  bending  toward  the  "common"  surface  atom , 
l.e.  bending  in  the  yz  plane.  The  results  are  shown  in  Table  9. 

The  coadsorption  effect  is  the  strongest  here ,  since  the  adsorbates 
are  very  close  to  each  other .  The  carbon  and  nitrogen  are  only 
1.90  A  apart.  As  in  the  case  of  high  coverage,  both  C-0  and  N-0 
have  larger  overlap  populations  than  when  they  are  alone  at  the 
surface.  The  C -surface  and  H -surface  interactions  are  weakened , 
with  respect  to  the  common  Rh  ( 2 )  atom ,  while  for  the  other  Rh 
atoms,  the  overlap  populations  are  slightly  larger  compared  with 
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Table  9.  The  coadsorption  oi  CO  and  NO  on  the  Rh(i00)  surface 


(see 

34)  with 

NO  bending  in  zz 

:  and  yz 

planes . 

xz 

plane 

yz 

plane 

* 

0° 

45° 

90° 

45° 

90° 

CO  Electron 

5 a 

1.67 

1.67 

1.67 

1.67 

1.65 

Densities 

2"x 

0.72 

0.72 

0.64 

0.61 

0.84 

2ny 

0.46 

0.41 

0.66 

0.46 

0.42 

NO  Electron 

50 

1.71 

1.60 

1.76 

1.61 

1.62 

Densities 

2wx(z) 

0.87 

1.05 

1  .01 

1.00 

1.06 

2iry(z) 

1.22 

0.98 

1  .  10 

0.96 

1.21 

Overlap 

c-o 

1.11 

0.96 

0.91 

0.99 

1.02 

Populations 

C-Rh  ( 2 ) 

0.47 

0.47 

0.36 

0.38 

0.14 

C-Rh( 1 ) 

0.63 

0.64 

0.50 

0.64 

0.65 

N-O 

0.96 

0.91 

0.93 

0.93 

0.84 

N-Rh  ( 2 ) 

0.41 

0.40 

0.33 

0.22 

-0.04 

H-Rh  ( 3 ) 

0.55 

0.47 

0.35 

0.48 

0.49 

C-H 

0.16 

0.04 

-0.05 

0.07 

0.01 

c-oHO 

-0.08 

0.05 

0.54 

0.04 

0.46 

0-0 

-0.04 

0.26 

-0.11 

0.25 

-0.10 

0H0-Rh(2) 

-0.07 

- 

- 

-0.10 

0.47 

Binding 
Energy  (eV) 

4.72 

0.79 

0.66 

-0.30 

1.59 

Pace  49 


the  single  adsorbate  studies.  Also,  there  is  a  certain  c-N 
attractive  interaction  which  occurs  primarily  through  their  n 
orbitals,  including  lw  orbitals.  The  carbon  and  nitric  oxide 
oxygen  do  not  interact,  since  they  are,  m  this  geometry, 
still  quite  far  away  from  each  other (2. 2 6  A). 

Let  us  now  discuss  the  HO  bending  in  the  xz  plane.  For 
#=45°,  the  C-o  and  H-0  overlap  populations  are  further  reduced. 
While  the  overall  C -surface  interaction  remains  essentially  the 

same ,  the  N  -surface  interaction  is  reduced .  The  weakening  of  the 
C-0  bond  is  mainly  due  to  the  oxygen -oxygen  interaction  which 
presumably  goes  through  the  lir  orbitals  (not  shown  in  the  Table 
9)  and  some  through  the  HO  o  orbitals.  We  still  have  some 
positive  C-H  overlap  populations  and  a  small  positive  carbon - 
oxygen  (from  HO)  overlap  population  of  0.05,  as  the  c-oHO  distance 
is  now  1.56  a.  When  we  bend  the  HO  further  #  =  90°,  this  distance 
is  reduced  to  1.36  A  which  results  in  a  high  C-Oj|0  interaction 

with  an  overlap  population  of  0.54.  The  C-0  overlap  population 

is  reduced  to  0.91.  The  H-0  overlap  population  (0.93)  is  now 
larger  than  for  #=45°,  due  to  the  smaller  H -surface  interaction, 

but  smaller  with  respect  to  HO  alone  at  the  surface  in  the  same 
geometry .  This  geometry  is  showing  a  strong  tendency  for  COg 
formation  and  some  small  effect  to  reduce  the  H-0  bond  strength, 
but  the  energetics  of  these  changes  do  not  really  favor  them . 

The  bending  of  HO  to  45°  in  the  yz  plane,  toward  the  surface 
Rh(S)  atom  is  energetically  unfavorable.  The  trends  m  overlap 
population  are  very  similar  to  those  of  bending  in  the  xz  plane, 
except  for  the  H  -surface  overlap  populations.  They  change 
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logically  and  the  H-Rh(2)  overlap  population  decreases 
significantly  and  N  -Rh  ( 3 )  decreases  also .  thougn  not  as  much  . 

Both  are  much  smaller  than  the  corresponding  values  when  HO  is 

alone  at  the  surface.  When  HO  is  lying  down  m  this 

arrangement ,  the  carbon  and  oxygen  from  HO  are  again  very  close  to 
each  other  (1.36  A ) ,  only  now  the  Rh  ( 2 )  surface  atom  is  "between" 
them .  Again  4o  and  in  CO  and  HO  orbitals  are  strongly  involved 
in  bonding.  There  is  a  depopulation  of  the  HO  lx  orbitals  and  a 
greater  donation  into  the  2x  •  s .  Hence ,  the  H  -0  overlap  population 
is  reduced .  The  C  -o  overlap  population  is  larger  than  when 

#=45°,  but  still  much  lower  than  when  CO  is  alone  at  the  surface. 

The  H-Rh  ( 2 )  overlap  population  is  negative  and  the  overlap 
population  between  carbon  and  Oho  strongly  increases.  The  CO- 
surface  interaction  is  further  reduced,  with  C-Rh(2)  overlap 
population  of  only  0.13,  so  that  the  CO  is  weaKly  bound  to  the 

surface  now.  If  we  have  a  situation  in  which  weaKly  bound  HO 

species  (lying  down  or  highly  inclined)  at  the  surface  approache 

coadsorbed  COs,  in  this  geometry,  then  COg  and  H  species  are 
possible  products ,  with  COg  being  weaKly  bound  to  the  surface . 

Finally ,  we  suggest  a  possible  geometry  where  both  CO  and  ho 
are  in  bridging  position  and  bound  to  adjacent  Rh  atoms  as  shown 
in  35.  The  bending  of  NO  is  in  the  xz  plane  here,  toward  the 

Rh ( 2 )  atom.  The  results  are  given  m  Table  10. 

The  coadsorption  effect  is  not  strong  (for  P~o°)  but 
notable.  The  occupation  of  the  CO  2ir*  orbitals  is  larger  than 

when  CO  is  alone  at  the  surface  (see  Table  l ) ,  but  the  CO  overlap 
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Table  10.  CO  and  HO  coadsorbed  on  the  Rh(iOO)  surface  in  the 


geometry 


CO  Electron 
Densities 


NO  Electron 
Densities 


Overlap 

Populations 


Binding 
Energy  (ev ) 


shown  in  35,  with  NO  bending  in  the  xz  plane. 


0° 

45° 

90° 

50 

1.65 

1.66 

1 .65 

2irx 

0.77 

0.74 

0.70 

2lfy 

0.44 

0.45 

0.42 

5o 

1.72 

1 .72 

1.44 

2wx(z) 

1.23 

0.73 

1.32 

2ffy 

0.61 

0.93 

0.96 

C-O 

1  .  10 

1  .  10 

1.06 

C-Rh  ( 2 ) 

0.62 

0.54 

0.26 

c-Rh(i ) 

0.62 

0.61 

0.62 

N-O 

0.94 

1.06 

0.85 

N  -Rh  ( 2  ) 

0.54 

0.31 

-0.02 

N-Rh  ( 3 ) 

0.54 

0.53 

0.45 

C-N 

0.01 

-0.02 

-0.03 

C“°N0 

0.00 

0.02 

0.33 

0-0 

0.00 

-0.06 

-0.10 

0H0-Rh(2) 

-0.06 

-0.09 

0.56 

4.96 

5.53 

4.62 

Page  52 


35 

population  is  somewhat  larger.  This  is  due  to  a  weaker  CO- 
surface  interaction.  The  same  is  true  for  NO  (see  Table  3). 
Again ,  this  is  consistent  with  the  experiments  which  suggest  that 
at  high  coverages ,  the  stretching  frequencies  shift  to  a  larger 
values  and  the  desorption  temperatures  decrease . 

For  the  bending  angle  of  45°,  both  CO-  and  NO-  surface 
interactions  are  reduced .  This  effect  is  greater  for  NO  and  Is 
reflected  In  a  lower  occupation  of  the  2s"  orbitals  and  larger 
N-o  overlap  population.  We  cannot  see  significant  interaction 
of  the  oxygen  from  NO  with  the  carbon  and  Rh(2)  atoms, even  though 
the  distances  of  2.04  and  2.22  A,  respectively ,  suggest  that  some 
interaction  might  occur.  When  NO  Is  lying  down,  these  distances 
are  reduced  to  1.54  and  1.36  A,  respectively .  The  CO-surface 
interaction  is  reduced  and  the  nitrogen  loses  its  interaction 
with  the  surface  Rh  ( 2 )  atom  .  The  N  -surf ace  interaction  here  is 
smaller  than  in  the  case  of  individual  NO  adsorption.  The  N-0 
overlap  population  is  strongly  reduced  and  is  smaller  than  when 
NO  is  alone  at  the  surface  in  the  same  geometry  (0.89 -Table  5). 
As  it  was  proposed  earlier,  this  may  be  a  candidate  for  the 
NO  dissociation  mechanism.  C-0  overlap  population  is  also 
further  reduced.  The  oxygen -oxygen  overlap  population  is 
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negative  and  oxygen  from  NO  is  now  more  weakly  bound  to  the 
Rh(2)  atom  than  in  the  case  of  individual  NO  adsorption 
(0.71  -Table  5 ) .  CO  can  interact  with  lying  down  NO  species  at 
the  surface  in  this  geometry ,  the  trends  for  COg  formation  and 
N-0  bond  breaking  are  evident. 

To  conclude  this  section  we  can  point  out  that  most  of 
the  NO  bending  mechanisms  do  not  favor  the  CO/NO  reaction , 
but  that  there  are  several  arrangements ,  with  NO  lying  down  at 
the  surface  which  would  lead  to  N-0  bond  braking  and  COg 
formation .  So ,  perhaps ,  we  do  not  need  the  NO  dissociation  as  a 
necessary  precondition  f o  this  reaction  to  occur .  In  that  sense 
we  could  support  the  conclusion  of  Mendershot  and  Hansen39  that 
the  NO  dissociation  is  not  important  under  their  reaction 
conditions.  On  the  other  hand  we  find  the  strongest  effects  if 
NO  is  lying  down  along  the  bridged  surface  Rh-Rh  bond  (as  in  34 
and  35),  i.e.  the  only  possible  candidate  for  NO  dissociation 
mechanism  in  the  frame  of  our  bending  scheme .  But ,  we  were  able 
to  produce  the  desired  trends,  though  to  a  smaller  extent,  in 
some  other  geometries ,  as  in  34  when  NO  is  lying  down  with  the 
oxygen  being  close  to  a  four -fold  hollow  site  (as  in  26),  or 
even  in  33  with  the  bending  angle  greater  than  90°. 

Coadsorption  of  CO  and  Oxygen 

Up  to  now  in  all  the  coadsorption  geometries ,  we  applied  the 
bending  model  for  NO  dissociation .  We  have  found  several 
geometries  of  lying -down  NO  with  which  the  CO  could  readily  react 
to  form  COg.  Now  we  would  like  to  briefly  investigate  the 
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coadsorption  of  CO  and  oxygen  as  the  interaction  of  such  species 
is  suggested  to  be  the  rate  determining  step  for  the  CO/NO 
reaction  on  Rh ( 1 1 1 ) .  The  adsorption  of  oxygen  on  Rh  surfaces 
has  been  investigated  experimentally2  • 55  • 61 .  On  the  Rh(iOO) 
surface"  oxygen  is  dissociatlvely  adsorbed  m  three  0  states, 
and  there  is  a  low  temperature  molecular  adsorption  state . 

Here  we  assumed  that  Og  has  dissociated  and  that  the  oxygen  are 
found  in  four -fold  hollow  sites.  The  Rh-0  distance  of  2.05  A  was 

taken  from  the  corundum  structure  of  RhgOj  system62. 

Keeping  the  oxygen  in  the  four -fold  hollow  we  have  two 

cu 

possible  ways  to  obtaln^close  C(0)-0  contact.  The  first  is  an 

on -top  sites  as  m  36,  and  the  second  is  with  CO  m  the  two -fold 
bridging  position,  37.  The  C-0(ads.)  distances  are  2.29  A  and 
1.53  A,  respectively  .The  results  for  these  sites  are  presented  in 
Table  1 1 . 


36  37 

Let  us  first  discuss  the  geometry  where  the  CO  is  on  -top ,  36 . 

When  coadsorbed  with  the  oxygen,  the  degeneracy  of  CO  »  orbitals 
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Table  11.  The  coadsorption  of  CO  and  o  on  the  Rh(iOO)  surface 
in  the  geometries  shown  in  36  and  37. 


36 

CO  or  0 
alonea 

37 

CO  or  0 
alonea 

CO  Electron 
Densities 

5o 

1.76 

1.78 

1.61 

1.72 

11Tx 

1.99 

1.99 

1.90 

1.96 

Iffy 

1.99 

1.99 

1.95 

1.95 

50 

1.64 

1.67 

1 .58 

1.69 

2™x 

0.32 

0.36 

0.61 

0.47 

2ffy 

0.32 

0.36 

0.73 

0.74 

O  electron 
densities 

s 

1.79 

1.78 

1.72 

1.78 

Px 

1.46 

1 .49 

1.30 

1.49 

Py 

1.46 

1  .49 

1.45 

1.49 

Pz 

1.92 

1.92 

1.75 

1.92 

Overlap 

Populations 

C-0 

1.26 

1.23 

1.04 

1  .  10 

C-Rh(l) 

0.82 

0.66 

0.44 

0.62 

Oa-Rh(  1 ) 

0.20 

0.23 

0.  12 

0.23 

Oa-Rh(2) 

0.23 

0.23 

0.24 

0.23 

c-oa 

0.01 

* 

0.47 

a  CO  in  on -top  or 

bridging  sitees ,  singly  adsorbed 

in  a  p  ( 2x  2  ) 

structure;  0  m 

four -fold 

hollow  site , 

singly  adsorbed  in  a 

p ( 2x2 )  structure. 
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is  removed ,  but  the  differences  between  these  orbitals  are  small . 


Also,  the  oxygen  px  and  py  orbitals  lose  their  degeneracy  and 
that  effect  is  slightly  greater.  CO  2i“  and  5c  orbitals  are 

somewhat  depopulated  when  compared  with  CO  being  alone  at  the 
surface .  The  same  is  true  for  oxygen  px  and  p7  atomic  orbitals . 

The  c-0  overlap  population  increases .  This  geometry  does  not 
provide  for  efficient  C-0  (ads.)  interaction  as  the  c-0  distance 
is  prohibitively  long .  Instead ,  we  see  repulsive  interaction  of 
coadsorbates  and  weaker  adsorbate  surface  bonding . 

In  the  case  where  CO  is  in  the  bridging  position,  the 
coadsorbates  are  closer  to  each  other ,  and  CO  tr  and  oxygen  p 

orbitals  are  collnear ,  so  it  is  logical  to  expect  more 
Interaction .  Indeed ,  the  changes  in  electron  densities  and 
overlap  populations  are  much  greater  than  in  the  previous  case. 
The  C-O  overlap  population  is  reduced  and  we  gain  a  positive 
C-0(ads.)  overlap  population  of  0.46.  At  the  same  time 
the  C-Rh(l)  overlap  population  is  reduced  as  are  the  overlap 
populations  between  the  adsorbed  oxygen  and  surface  rhodium.  The 
geometry  looks  promising  for  the  formation  of  CO  2  molecules  at  the 
surface  and  subsequent  desorption  as  well . 
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Appendix 

In  this  paper,  all  performed  computations  are  tight -binding 
band  calculations  of  the  extended  Huckel  type.  The  H^'s  for 
rhodium  were  obtained  by  charge  iteration  on  the  bulk  metal 
(with  the  experimental  fee  geometry)  using  Grays  equation63, 
with  A ,  B  and  C  parameters  taken  from  Ref .  64 .  The  carbon 
Hi±>s  were  taken  from  an  earlier  work  of  this  group12  and 
were  in  good  agreement  with  iterated  values .  The  oxygen 
Hu’S  were  obtained  by  charge  iteration  on  CO  adsorbed  on 
both  sides  of  the  rhodium  three -layer  slab,  keeping  the  c 
and  Rh  Hu’S  constant,  with  the  low  coverage  of  0.25  on  both 
sides  to  avoid  the  adsorbate  interactions.  Hltrogen  HiA >s 
were  obtained  for  HO  adsorbed  in  the  same  way  as  CO,  with  the 
same  procedure,  with  oxygen  and  rhodium  parameters  being 
constant.  Oxygen  parameters  were  previously  determined 
parameters  from  CO  adsorption .  B  and  C  iteration  parameters 
for  C,  H  and  o  atoms  were  taken  from  Ref.  65.  Extended  Huckel 
parameters  for  all  atoms  used ,  are  listed  in  Table  i  2 . 

All  the  results  listed  in  this  paper  are  from  the  single  - 
face  adsorption  calculations ,  using  the  three  -layer  slab  model . 
The  geometrical  parameters  for  coadsorption  are:  C-0,  l .  15;  H-o, 
1.15;  Rh-Rh,  2.6693;  Rh-C,  1.90;  Rh-H,  1.90;  and  Rh-O,  2.05 
(all  values  in  A).  The  k -point  sets  varied  from  49  to  12 
points  in  order  to  keep  the  total  number  of  states  for  each 
geometry  similar. 
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Table  1 2 .  Extended  HucKel  parameters 


Orbital 

HU  <eV) 

5* 

cla 

Rh  5s 

-7.31 

2.13 

RH  5p 

-3.39 

2.10 

Rh  4d 

-10.35 

4.29 

1.97 

0.5607 

C  2s 

-16.20 

1.63 

C  2p 

-9.50 

1.63 

N  2S 

-23.95 

1.95 

N  2p 

-10.95 

1.95 

O  2s 

-27.61 

2.26 

O  2p 

-11.01 

2.26 

a  Contraction  parameters  used  in  tne  double-  j  expansion. 


c2a 


0.5665 
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